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Abstract 


This  report  describes  a  computational  study  undertaken  to  consider  the 
aerodynamic  effect  of  small  tiny  jets  as  a  means  to  provide  the  control  authority 
needed  to  maneuver  a  projectile  at  low  subsonic  speeds.  Scalable  Navier-Stokes 
computational  techniques  have  been  used  to  obtain  numerical  solutions  for  the 
jet-interaction  flow  field  for  a  projectile  at  subsonic  speeds.  Computed  results 
have  been  obtained  at  low  subsonic  speeds  at  0°  and  4°  angle  of  attack.  Both 
steady  and  unsteady  jets  have  been  considered.  For  comparison  purposes,  a 
jet-off  case  was  also  computed.  Qualitative  flow  field  features  show  the 
interaction  of  jets  with  the  free  stream  flow.  Numerical  results  show  the  effect  of 
the  jet  locations  and  sizes  on  the  flow  field  and  surface  pressures,  and  hence  on 
title  aerodynamic  coefficients.  Unsteady  jet  results  have  been  obtained  for  a 
two-dimensional  (2-D)  jet  flow  and  compared  with  experimental  data  for 
validation.  Some  results  obtained  with  an  unsteady  jet  for  the  subsonic  projectile 
are  included.  These  numerical  results  are  being  assessed  to  determine  if  small 
tiny  jets  can  be  used  to  provide  the  control  authority  needed  for  maneuvering 
munitions  in  lieu  of  canards  and  fins. 
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1.  Introduction 


The  prediction  of  aerodynamic  coefficients  for  projectile  configurations  is 
essential  in  assessing  the  performance  of  new  designs.  Accurate  determination 
of  aerodynamics  is  critical  to  the  low-cost  development  of  new  advanced  guided 
projectiles,  rockets,  missiles,  and  smart  munitions.  Fins,  canards,  and  jets  can  be 
used  to  provide  control  for  maneuvering  projectiles  and  missiles.  The  flow  fields 
associated  with  these  control  mechanisms  for  U.S.  Army  weapons  are  complex 
involving  three-dimensional  (3-D)  shock-boundary  layer  interactions, 
jet-interaction  with  the  free  stream  flow,  and  highly  viscous-dominated 
separated  flow  regions.  The  jet  interference  extends  over  significant  portions  of 
the  projectile  and  must  be  modeled  correctly.  For  missiles,  jet  thrusters  have 
been  studied  over  a  number  of  years  to  provide  high-speed  aerodynamic  control. 
These  thrusters  interact  with  the  surrounding  flow  field;  again,  the  resulting  jet 
interaction  flow  field  is  complex.  Recently,  several  studies  have  shown  that 
small  tiny  microelectromechanical  system  (MEMS)  jets  can  significantly  alter  the 
flow  field  and  pressure  distributions  for  airfoils  and  cylinders.  The  present 
analysis  involves  the  use  of  these  tiny  jets  for  projectile  aerodynamic  control. 
The  emphasis  in  the  present  research  is  to  provide  insight  into  the  interaction  of 
these  jets  with  the  free  stream  flow  and  to  determine  the  feasibility  of  these  jets 
for  aerodynamic  control  of  a  subsonic  projectile.  Both  computational  and 
experimental  data  for  these  jet  interactions  are  very  limited.  Simple  theories 
cannot  predict  the  complex  flow  fields  associated  with  die  jet  interaction  and 
experimental  tests  are  very  expensive.  To  help  reduce  experimental  costs, 
computational  fluid  dynamics  (CFD)  is  being  used  to  predict  these  complex 
flows  and  to  provide  detailed  pressure  and  force  and  moment  data.  The 
advanced  CFD  capability  used  here  solves  the  Navier-Stokes  equations  and 
incorporates  the  Chimera  technique.  This  numerical  capability  has  been  used 
successfully  to  determine  the  aerodynamics  on  a  number  of  complex  military 
systems  [1-3]  of  interest  to  the  U.S.  Army. 

This  report  describes  the  application  of  CFD  to  a  40-mm  grenade  with  MEMS 
jets.  Figure  1  shows  a  model  of  this  projectile.  The  problem  involves  3-D  flow 
computations  on  a  medium-caliber  subsonic  grenade  with  multiple  jet  holes. 
However,  the  scope  of  this  report  is  limited  to  a  single  jet  hole.  By  varying  the 
size  of  the  jet  hole  as  well  as  its  location  on  the  projectile  body,  a  range  of 
computational  results  has  been  provided.  Numerical  flow  field  computations 
have  been  made  for  both  steady  and  unsteady  jets  at  a  low  subsonic  speed.  The 
unsteady  jet  case  for  this  subsonic  projectile  includes  modeling  the  inside  of  the 
jet  cavity.  The  purpose  is  to  determine  the  resulting  aerodynamic  forces  on  the 
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Figure  1.  Model  of  the  40-mm  grenade. 

grenade  when  the  jet  flow  is  released  from  various  size  cavities,  at  each  of  the 
two  locations  (one  on  the  nose  of  the  projectile  and  one  on  the  cylinder). 


2.  Solution  Technique 


The  complete  set  of  time-dependent,  Reynolds-averaged,  thin-layer 
Navier-Stokes  equations  is  solved  numerically  to  obtain  a  solution  to  this 
problem.  The  numerical  technique  used  is  an  implicit,  finite  difference  scheme. 
Time-accurate  calculations  are  made  to  numerically  simulate  a  medium-caliber 
grenade  with  a  jet  hole  in  various  locations  on  the  surface  of  the  projectile. 

2,1  Governing  Equations 

The  complete  set  of  3-D,  time-dependent,  generalized-geometry, 
Reynolds-averaged,  thin-layer  Navier-Stokes  equations  is  solved  numerically  to 
obtain  a  solution  to  this  problem  and  can  be  written  in  general  spatial 
coordinates  §,  q,  and  ^  as  follows  [4]: 

d Tq  +  d gF +  d nG  +  d  d ^S,  (1) 

where 
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£  =  $(x,  y,  z,  t) —longitudinal  coordinate, 

T]  =  r|(x,  y,  z,  t)  —  circumferential  coordinate, 

£  =  £(x/  Y'  z' t) —nearly  normal  coordinate,  and 
t  =  t— time. 


In  equation  1,  q  contains  the  dependent  variables  (density,  three  velocity 

components,  and  energy)  and  F  ,G ,  and  H  are  flux  vectors.  The  thin-layer 
approximation  is  used  here,  and  the  viscous  terms  involving  velocity  gradients 
in  both  the  longitudinal  and  circumferential  directions  are  neglected.  The 
viscous  terms  are  retained  in  the  normal  direction,  Z,,  and  are  collected  into  the 

vector  S .  In  the  wake  or  the  base  region,  similar  viscous  terms  [1]  are  also  added 
in  the  streamwise  direction,  An  implicit,  approximately  factored  scheme  is 
used  to  solve  these  equations.  For  computation  of  turbulent  flows,  the  turbulent 
contributions  are  supplied  through  an  algebraic  eddy  viscosity  turbulence  model 
developed  by  Baldwin  and  Lomax  [5]  or  a  pointwise  turbulence  model  [6]. 

2.2  Numerical  Algorithm 

The  implicit,  approximately  factored  scheme  for  the  thin-layer  Navier-Stokes 
equations  using  central  differencing  in  the  q  and  £  directions  and  upwinding  in  t; 
is  written  in  the  following  form  [7]: 

[l  +  ibhSb4  (A+y+ibhS<Cn-ibh  Re"1  gc  J-'  MnJ-ibD, ,  ^  J 


x[  I  +  ibhS/  ( A-y+ibhSnBn~ibDi  |,  ]  A Q” 


=ibAt  {  <?*>[  (f+  y-F- 


]  +  S/[(F-)”-F_l]  +  Sr,(Gn-Gx) 


(2) 


+  Se  (Hn-Hx)-Re-'6((Sn-S„)  }-ibDe(Q” -QJ  , 


where  h  =  At  or  (At)/ 2  and  the  free-stream  base  solution  is  used.  The  free-stream 
fluxes  are  subtracted  from  the  governing  equation  to  reduce  the  possibility  of 
error  from  the  free-stream  solution  corrupting  the  converged  solution.  Here,  8  is 
typically  a  three-point  second-order  accurate  central  difference  operator,  6  is  a 
midpoint  operator  used  with  the  viscous  terms,  and  the  operators  Sb  and  S$ 
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are  backward  and  forward  three-point  difference  operators.  The  flux  F  has 
been  eigensplit,  and  the  matrices  A,B,C,  and  M  result  from  local  linearization 
of  the  fluxes  about  the  previous  time  level.  Here,  J  denotes  the  Jacobian  of  the 
coordinate  transformation.  Dissipation  operators  De  and  D,  are  used  in  the 
central  space  differencing  directions.  The  smoothing  terms  used  in  the  present 
study  are  of  the  form 


and 


where 


£>,|,  =  (A0/-' 

s2  8p(B) (5 S  +e4S  p(B)  s  3 

1  +  P 

[s2Sp(B)j3S+ 2.5  s45p(B)5  ]| 

fi-  , 

1  (  1  +  £2  )  P  | 

J, 


(3) 


(4) 


and  where  p(B)  is  the  true  spectral  radius  of  B.  The  idea  here  is  that  the  fourth 
difference  will  be  tuned  down  near  shocks  (e.g.,  as  j?  gets  large,  the  weight  on  the 
fourth  difference  drops  down  while  the  second  difference  tunes  up). 

2.3  Boundary  Conditions 

For  simplicity,  most  of  the  boundary  conditions  have  been  imposed  explicitly. 
The  no-slip  boundary  condition  is  used  on  the  grenade  surface;  it  has  been 
modified,  however,  to  allow  a  jet  flow  of  variable  size  and  location.  The  pressure 
at  the  wall  is  calculated  by  solving  a  combined  momentum  equation.  A 
symmetry  boundary  condition  is  imposed  at  the  circumferential  edges  of  the 
grid.  At  the  centerline  axes,  a  polar  boundary  condition  is  used  at  the  nose  of  the 
projectile,  while  a  collapsed  axis  boundary  condition  is  used  at  the  centerline  in 
back  of  the  projectile.  Boundary  conditions  are  not  applied  at  the  outer 
boundary,  which  is  far  enough  from  the  body  to  allow  for  free  stream  flow. 


3.  Chimera  Composite  Grid  Scheme 


The  chimera  overset  grid  technique  [8-10]  involves  generating  independent  grids 
about  each  component  and  then  oversetting  them  onto  a  base  grid  to  form  the 
complete  model.  This  procedure  reduces  a  complex  single  or  multi-body 
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problem  into  a  number  of  simpler  subproblems.  An  advantage  of  the  overset 
grid  technique  is  that  it  allows  computational  grids  to  be  obtained  for  each  body 
component  separately  and  thus  makes  the  grid  generation  process  easier. 
Because  each  component  grid  is  generated  independently,  portions  of  one  grid 
may  lie  within  a  solid  boundary  contained  within  another  grid.  Such  points  lie 
outside  the  computational  domain  and  are  excluded  from  the  solution  process. 
Equation  2  has  been  modified  for  chimera  overset  grids  by  the  introduction  of 
the  flag  ib  to  achieve  just  that.  This  if,  array  accommodates  the  possibility  of 
having  arbitrary  holes  in  the  grid.  The  if,  array  is  defined  so  that  if,  =  1  at  normal 
grid  points  and  if,  =  0  at  hole  points.  Thus,  when  if,  =  1,  equation  2  becomes  the 
standard  scheme.  The  set  of  grid  points  that  forms  the  border  between  the  hole 
points  and  the  normal  field  points  is  called  inter-grid  boundary  points.  These 
points  are  updated  by  interpolating  the  solution  from  the  overset  grid  that 
created  the  hole.  Values  of  the  if,  array  and  the  interpolation  coefficients  needed 
for  this  update  are  provided  by  a  separate  algorithm  [10]. 


4.  Model  Geometry  and  Computational  Grid 


The  computational  model  consists  of  a  40-mm  grenade,  1.76  calibers  in  length. 
The  mesh  for  this  model  consists  of  a  single  grid,  containing  approximately 
750,000  grid  points.  Figure  2  shows  a  cross-sectional  view  of  the  longitudinal 
grid  system  for  this  configuration.  The  grid  dimensions  are  151  x  61  x  80  in  the 
longitudinal,  circumferential,  and  normal  directions  respectively.  Initially,  the  jet 
holes  were  modeled  by  modifying  the  boundary  conditions  at  specific  indices  on 
the  surface  of  the  grenade.  The  jet  width  is  1.12  mm  (0.0281  calibers).  In  the 
circumferential  direction,  the  jet  was  modeled  for  various  size  openings: 
6°  (small),  45°  (medium),  and  90°  (large).  In  Figures  3  and  4,  the  jet  openings  are 
shaded  in  black.  Figure  3  shows  their  location  on  the  nose  of  the  grenade. 
Similar  jet  holes  were  also  modeled  on  the  cylinder  of  the  grenade  at 
X/D  =  4.6  mm  (1.49  calibers),  as  shown  in  Figure  4.  This  mesh  was  used  for  the 
steady  jet  computations.  For  the  jet-off  cases,  a  similar  axisymmetric  grid  was 
used,  with  only  three  planes  in  the  circumferential  direction. 

For  the  unsteady  jet  calculations,  the  computational  mesh  for  the  medium  jet  was 
modified  to  include  the  jet  cavity  inside  the  projectile.  The  depth  of  the  cavity  is 
10  mm  (0.25  calibers).  The  original  mesh  used  for  the  projectile  body  was  refined 
by  increasing  the  number  of  points  in  the  axial  direction  to  187  and  by  clustering 
the  points  in  the  areas  of  the  nose  and  cylinder  jets.  In  addition,  a  small  mesh 
was  generated  to  model  the  jet  cavity  inside  die  projectile  body.  This  mesh 
consists  of  approximately  7,000  grid  points:  11  x  16  x  41  in  the  axial, 
circumferential,  and  normal  directions.  Consequently,  the  total  number  of  grid 
points  increased  to  approximately  920,000.  Figure  5  shows  a  computational 
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Figure  2.  Computational  model  and  grid  system. 


Figure  3.  Location  of  nose  jets  on  the  computational  model. 


I 

small  medium  large 

Figure  4.  Location  of  cylinder  jets  on  the  computational  model. 


Figure  5.  Computational  model  showing  jet  cavity. 

model  of  the  projectile  with  the  jet  cavity.  Since  the  Chimera  technique  was  to  be 
used  for  the  unsteady  jet  computations,  the  jet  cavity  mesh  could  be  moved  to 
various  locations  if  necessary.  PEGSUS  [11],  a  mesh  interpolation  code,  was  used 
to  provide  the  necessary  intergrid  information. 
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5.  Results 


CFD  computations  have  provided  interesting  information  about  the  40-mm 
projectile  flow  field.  The  results  provide  insight  to  the  problem  and  demonstrate 
the  capability  of  the  Zonal  Navier-Stokes  flow  solver  ZNSFLOW  [12]  to  provide 
flow  field  solutions  for  a  projectile  with  jets. 

5.1  Jet-Off 

The  preliminary  effort  was  to  compute  axisymmetric  solutions  at  Mach  numbers 
0.15,  0.2,  and  0.25,  at  0°  angle  of  attack  and  wall  temperature  of  530°  Rankin. 
These  calculations  required  only  small  amounts  of  time  and  memory  and  were 
completed  on  an  SGI  Onyx  system  at  the  U.S.  Army  Research  Laboratory  (ARL) 
Major  Shared  Resource  Center  (MSRC)  site  at  Aberdeen  Proving  Ground,  MD. 
Figure  6  shows  the  convergence  history  of  the  drag  force  coefficient  (CDO)  for 
the  jet-off  cases.  Computed  pressure  contours  for  these  cases  are  shown  in 
Figure  7  and  show  the  expected  features  as  the  Mach  number  increases  from  0.15 
to  0.25  (top  to  bottom).  Figure  8  shows  that  the  change  in  drag  coefficient  as  a 
function  of  Mach  number  is  rather  small  in  the  range  of  Mach  numbers 
considered.  As  expected,  the  pressure  (wave)  drag  (CDP)  is  the  largest 
contributor;  the  skin  friction  drag  (CDV)  component  is  small.  Velocity  vectors  in 
the  base  region,  as  seen  in  Figure  9,  show  the  flow  to  be  similar  for  the 
computations.  Consequently,  a  Mach  number  of  0.25  was  selected  for  the  jet-on 
computations. 


Figure  6.  Convergence  history,  jet-off. 
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Figure  7.  Pressure  contours,  jet-off;  Mach  numbers  0.15, 0.2,  and  0.25  (top  to  bottom) 
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Figure  8.  Mach  number  vs.  drag  coefficients,  jet-off. 


Figure  9.  Velocity  vectors,  base  region,  jet-off. 
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5.2  Steady  Jet 

For  this  study,  steady-state  numerical  computations  were  performed  with  the  jet 
on,  varying  the  size  and  location  of  the  jet  holes  on  a  medium  size  grenade.  A  jet 
velocity  of  Mach  0.125  (one-half  the  free  stream  velocity)  was  used. 
Computations  were  completed  for  a  Mach  number  of  0.25,  and  at  0°  and  4° 
angles  of  attack.  Sea-level  atmospheric  flight  conditions  were  used.  Each  3-D 
calculation  required  approximately  6.6  million  words  of  memory,  and  each  case 
used  an  average  of  25  hr  of  computer  time  on  the  Silicon  Graphics  Origin  2000 
suite  of  supercomputers  at  the  ARL  MSRC.  Resources  used  for  the  axisymmetric 
cases  were  minimal. 

Axisymmetric  calculations  were  done  at  0°  angle  of  attack  for  the  nose  jet  and  the 
cylinder  jet.  Computational  results  for  the  jet-on  cases  were  compared  with  that 
of  the  jet-off  case.  Figure  10  shows  that  the  nose  jet  slightly  affects  the  surface 
pressure  distribution  on  the  nose  section  of  the  projectile,  while  the  cylinder  jet 
has  a  stronger  effect  both  upstream  and  downstream  of  the  jet  location.  The  drag 
coefficients  in  Figure  11  show  a  large  decrease  due  to  the  nose  jet,  while  the 
cylinder  jet  increases  drag  only  slightly.  The  pressure  contours  shown  in  Figure 
12  reveal  only  slight  changes  in  the  flow  field  at  the  nose  of  the  projectile.  The 
changes  are  much  more  dramatic  for  the  cylinder  jet  and  are  shown  in  Figure  13. 
The  differences  in  the  flow  field  upstream  as  well  as  downstream  of  the  jet  are 
clearly  visible,  especially  for  the  large  jet  case. 

Computed  results  have  been  obtained  for  3-D  steady  jet-on  cases  where  a  jet  was 
placed  either  on  the  nose  or  the  cylinder  of  the  projectile.  By  specifying  indices 
on  the  projectile  surface,  jet  flow  boundary  conditions  were  provided  for  jets  of 
three  sizes  (6°,  45°,  and  90°)  on  the  nose  of  the  grenade  as  well  as  on  the 
cylindrical  section  of  the  projectile  (see  Figures  3  and  4).  A  series  of 
computations  for  the  various  size  jets  and  locations  was  performed  for  both 
0°  and  4°  angle  of  attack.  Surface  pressure  data  was  extracted  at  the  0°  plane  for 
all  cases.  Figure  14  shows  a  comparison  of  surface  pressure  coefficient  for  all 
three  nose  jets.  Again  the  differences  are  mainly  on  the  nose  of  the  projectile. 
For  the  cylinder  jets,  the  differences  are  quite  dramatic,  both  upstream  and 
downstream  of  the  jet,  as  seen  in  Figure  15.  A  top-side  view  of  surface  pressure 
contours  on  the  body  are  shown  in  Figure  16  for  the  jet-off  case  and  all  three 
cylinder  jets.  (The  jet  location  is  shaded  in  gray.)  The  effects  of  the  jet  are  seen 
quite  clearly.  Longitudinal  pressure  contours  of  the  flow  fields  shown  in  Figure 
17  also  demonstrate  this  effect. 

Force  and  moment  coefficients  were  computed  for  all  solutions.  CN  is  normal 
force,  CA  is  axial  force  (drag),  and  CMP  is  the  pitching  moment.  Table  1  contains 
force  and  moment  data  for  the  nose  jet  locations  and  shows  that  there  are  only 
minimal  differences  for  the  various  size  jets.  Table  2  shows  a  slightly  stronger 
difference  for  the  cylinder  jet  locations,  especially  for  the  large  jet. 
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Figure  10.  Surface  pressure  coefficient  vs.  axial  position. 
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Figure  11.  Drag  coefficient  vs.  axial  position. 
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Figure  13.  Pressure  contours  in  cylinder  area,  (a)  jet-off  and  (b)  jet-on. 

5.3  Unsteady  Jet 

For  computational  validation  of  the  unsteady  jet  CFD  modeling,  an  isolated  two- 
dimensional  (2-D)  case  was  first  selected.  A  schematic  diagram  and  a  flow 
picture  obtained  from  the  experiment  for  this  isolated  2-D  jet  are  shown  in 
Figure  18.  Here  the  jet  width  is  0.5  mm  and  the  peak  jet  velocity  is  20  m/s.  In 
the  experiment,  the  synthetic  jet  was  formed  in  air  at  an  orifice  measuring 
0.5  mm  (width)  x  75  mm  (length).  The  jet  actuator  operates  at  a  frequency  of 
1,000  Hz  [13].  The  jet  is  synthesized  by  the  time-harmonic  motion  of  a  flexible 
diaphragm  in  a  sealed  cavity.  This  motion  results  in  both  positive  and  negative 
velocities  at  the  jet  exit  while  the  net  mass  flux  out  of  the  cavity  is  zero  during 
each  cycle.  In  the  computations,  unsteady  jet  boundary  conditions  were  applied 
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NOSE  JETS 


Figure  14.  Surface  pressure  coefficient— nose  jets. 


CYLINDER  JETS 


Figure  15.  Surface  pressure  coefficient— cylinder  jets. 

with  a  sinusoidal  variation  in  the  jet  velocity  with  a  peak  amplitude  of  20  m/s. 
Computed  velocity  and  vorticity  contours  are  shown  in  Figure  19.  The 
time-averaged  jet  centerline  velocity  over  many  cycles  of  unsteady  jet  CFD 
computations  is  compared  with  available  time-averaged  experimental  data  in 
Figure  20  and  is  found  to  be  in  reasonable  agreement.  As  shown  both  in  the 
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Medium  Jet 


Large  Jet 


Table  1.  Force  and  moment  coefficients  for  nose  jets,  a  =  4. 


Figure  18.  Schematic  and  flow  picture  of  a  2-D  jet  experiment. 


Figure  19.  Velocity  components  u,  v,  and  vorticity  for  an  unsteady  jet. 


experiment  and  the  computations,  the  time-averaged  centerline  velocity  is 
decreased  with  increasing  distance  away  from  the  jet  exit.  The  next  step  was  to 
apply  similar  unsteady  jet  boundary  conditions  to  the  subsonic  projectile  case. 
For  this  case,  the  jet  cavity  inside  the  projectile  was  included  (see  Figure  5).  The 
unsteady  jet  boundary  conditions  were  applied  at  the  bottom  wall  of  the  jet 
cavity.  Numerical  results  were  obtained  for  this  unsteady  jet  case  at  Mach  =  0.25 
and  zero  angle  of  attack.  The  jet  width  was  2.24  mm  and  the  peak  jet  velocity 
used  was  43  m/s  operating  at  a  frequency  of  1,000  Hz.  Computed  surface 
pressures  obtained  at  a  given  instant  in  time  are  shown  in  Figure  21.  It  also 
shows  particle  traces  coming  out  of  the  cavity.  The  pressure  field  both  upstream 
and  downstream  of  the  jet  is  affected  by  the  jet  flow  depending  on  whether  the 
flow  is  into  or  out  of  the  cavity.  Additional  qualitative  features  of  the  flow  field 
are  shown  in  Figure  22.  This  figure  shows  the  snapshots  of  the  computed 
velocity  vectors  in  the  vicinity  of  the  jet  exit  at  two  different  instants  in  time. 
Figure  22(a)  dearly  shows  the  external  free  stream  flow  going  into  the  jet  cavity 
whereas  Figure  22(b)  shows  the  flow  coming  out  of  the  cavity  and  interacting 
with  the  free  stream  flow.  Figure  22(b)  also  shows  a  region  of  separated  flow  just 
downstream  of  the  jet  exit.  The  resulting  surface  pressures  are  integrated  to 
obtain  the  aerodynamic  forces  and  moments.  The  computed  axial  force,  the 
normal  force,  and  the  pitching  moment  coefficients  are  shown  in  Figure  23  as  a 
function  of  time  for  a  cycle.  These  results  dearly  indicate  the  unsteady  nature  of 
the  flow.  Similar  results  were  obtained  for  this  projectile  at  an  angle  of  attack  of 
4°  (see  Figure  24).  These  results  are  being  analyzed  to  determine  the  feasibility  of 
these  jets  to  provide  control  authority.  It  is  anticipated  that  multiple  jets  may  be 
required  to  provide  the  control  authority  needed  for  maneuvering  a  subsonic 
projectile. 


Figure  21.  Computed  surface  pressures,  M  =  0.25,  a  =  0,  unsteady  jet. 
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Force  and  Moment  Coefficients 


(a)  (b) 

Figure  22.  Velocity  vectors  at  two  instants  in  time  during  the  cycle,  M  =  0.25,  a  =  0, 
unsteady  jet. 
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Figure  23.  Force  and  moment  coefficients,  M  =  0.25,  a  =  0,  unsteady  jet. 
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Figure  24.  Force  and  moment  coefficients,  M  =  0.25,  a  =  0,  unsteady  jet. 


6.  Conclusion 


A  computational  study  has  been  undertaken  to  consider  the  aerodynamic  effect 
of  small  tiny  jets  as  a  means  to  provide  the  control  authority  needed  to  maneuver 
a  projectile  at  low  subsonic  speeds.  Computed  results  have  been  obtained  at 
subsonic  speeds  and  at  0°  and  4°  angle  of  attack.  Both  steady  and  unsteady  jets 
have  been  simulated.  Qualitative  flow  field  features  show  the  interaction  of  tiny 
jets  and  the  extent  of  their  influence  both  upstream  and  downstream  of  the  jet. 
The  CFD  results  also  show  the  effect  of  the  jet  locations  and  sizes  on  the  surface 
pressure  distribution.  The  unsteady  jet  results  obtained  for  a  2-D  jet  are 
compared  with  the  experimental  data  and  are  found  to  be  in  reasonable 
agreement.  The  unsteady  jet  has  been  applied  to  the  subsonic  projectile  and  is 
shown  to  have  some  effect  on  forces  and  moments  even  at  zero  degree  angle  of 
attack.  The  results  show  the  potential  of  CFD  to  provide  insight  into  the  jet 
interaction  flow  fields  and  provide  guidance  as  to  the  locations  and  sizes  of  the 
jets  to  generate  the  maximum  control  authority  for  maneuvering  smart 
munitions. 
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